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QCD and a Holographic Model of Hadrons 
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We propose a five-dimensional frameworlc for modeling low-energy properties of QCD. In the 
simplest three parameter model we compute masses, decay rates and couplings of the lightest mesons. 
The model fits experimental data to within 10%. The framework is a holographic version of the QCD 
sum rules, motivated by the anti-de Sitter/conformal field theory (AdS/CFT) correspondence. The 
model naturally incorporates properties of QCD dictated by chiral symmetry, which we demonstrate 
by deriving the Gell-Mann-Oakes-Renner relationship for the pion mass. 



Introduction. — QCD has eluded an analytic solution, 
despite extensive efforts applied to this problem in the 
past 30 years. Recently, the gravity/gauge, or anti- 
de Sitter/conformal field theory (AdS/CFT) correspon- 
dence Ul has revived the hope that QCD can be reformu- 
lated as a solvable string theory. So far, theories which 
can be solved using AdS/CFT techniques differ substan- 
tially from QCD, most notably by the strong coupling 
in the ultraviolet (UV) regime and the lack of asymp- 
totic freedom. Nevertheless, certain important proper- 
ties of QCD, such as confinement and chiral symmetry 
breaking, are present in many of these theories, and the 
gravity/gauge duality provides a new approach to study- 
ing the resulting dynamics. An important development 
in the prototypical example of A/" = 4 super Yang-Mills 
(SYM) theory has been the introduction of fundamental 
quarks using probe D7 branes The mesons that ap- 
pear in these theories behave in many ways similarly to 
the mesons in QCD 0,0. 

Inspired by the gravity/gauge duality we propose the 
following complementary approach. Rather than deform 
the SYM theory to obtain QCD 0, we start from QCD 
and attempt to construct its five-dimensional (5D) holo- 
graphic dual. In this Letter, we present an exploratory 
study of a simple holographic model of QCD. The field 
content of the 5D theory is chosen to reproduce holo- 
graphically the dynamics of chiral symmetry breaking in 
QCD, the boundary theory. The model has four free pa- 
rameters, one of which is fixed by the number of colors; 
the remaining three parameters can be fitted using three 
well-measured observables, e.g., the p meson mass, the 
pion mass, and the pion decay constant. The model then 
predicts other low-energy hadronic observables with sur- 
prisingly good accuracy. 

Such an approach is similar in spirit to the construc- 
tion of the QCD moose theory in Ref. where the holo- 
graphic description arises in the continuum limit of in- 
finitely many hidden local symmetries (see also Ref. 0|)- 
As in Ref. ja], vector meson dominance and QCD sum 
rules are natural consequences of our model. Hence, the 
success of the model is not coincidental, but a result of 



linking several proven approaches through the AdS/CFT 
correspondence. We expect the success of our model to 
diminish above roughly the scale given by the mass of the 
lightest isospin-carrying spin- 2 resonance, namely the 02 
(1318 MeV In particular, we are completely neglect- 
ing stringy physics which becomes important at higher 
energies, and we have not included in our description any 
modes with spin larger than 1. At this stage, we also ne- 
glect running of the QCD coupling, which is likely a poor 
approximation for a larger range of energies. While our 
model is too simple to provide a complete dual descrip- 
tion of QCD, its success seems to suggest that there is 
a quantitatively useful reformulation of QCD as a string 
theory in a higher-dimensional curved space. 

Field content. — Table HI illustrates the field content of 
our model. The choice of the 5D fields is dictated by a 
principle of the AdS/CFT correspondence: each opera- 
tor 0{x) in the 4D field theory corresponds to a field 
4>{x, z) in the 5D bulk theory. The 5D theory dual 
to QCD should, therefore, contain an infinite number 
of fields corresponding to the infinite number of oper- 
ators in QCD. There is, however, a small number of 
operators that arc important in the chiral dynamics: 
the left- and right-handed currents corresponding to the 
S\]{N f)Lx'SU{N f)fi chiral fiavor symmetry, and the chi- 
ral order parameter (see Table P) . We shall include in 
our model only the 5D fields which correspond to these 
operators and neglect all other fields. 

The 5D masses ms of the fields ^2^, A°^^, and X are 
determined via the relation 9, 10] (A — p)(A -f p — 4) = 
TO5, where A is the dimension of the corresponding p- 
form operator — see Tabled We assumed here that these 
operators keep their canonical dimensions, which is true 

TABLE I: Operators/fields of the model 



4D: 0(x) 


5D: (f){x, z) 
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Aa 
ARii 
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3 
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qRqi 







3 


-3 
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only for the conserved currents. However, for the field 
X we could easily incorporate corrections to its classical 
dimension. The factor 1/z in Table ^ is dictated by the 
dimension of the operator qq, while the factor of 2 is of no 
physical significance and is chosen for later convenience. 

We shall choose the simplest possible metric for our 
model, namely, a slice of the anti-de Sitter (AdS) metric. 



ds^ — —{—dz^ + dx^dxi_i), < z < z„ 



(1) 



The fifth coordinate z corresponds to the energy scale, as 
higher energy (or momentum transfer Q^) QCD physics 
is reflected by the behavior of the fields closer to the 
AdS boundary z = 0: Q ^ 1/z. By virtue of the con- 
formal isometry of the AdS space, in such a model the 
running of the QCD gauge coupling is neglected in a win- 
dow of scales until an infrared (IR) scale Qm ^ ^/zm- To 
make the theory confining, one introduces an IR cutoff 
in the metric aX z = Zm where spacetime ends, in anal- 
ogy with the case of the cascading gauge theory studied 
in Ref. 0|. We shall call z — Zm the "infrared brane" 
and impose certain boundary conditions on the fields at 
z = Zm- Certainly, this is only a crude model of con- 
finement. Indeed, our model requires two dimensionful 
parameters related to chiral symmetry breaking, whereas 
in QCD there is only one. In addition, an UV cutoff can 
be provided by setting the boundary to z — e instead of 
z = Below we shall frequently use such a cutoff as a 
mathematical tool, but we shall always imply the limit 
of e — > for simplicity. 

5D action and chiral symmetry breaking. — The action 
of the theory in the bulk is 

^^{Fl + Ff,)] (2) 

where D^,X = d^,X - iAL^^X + iXAa^, Alm = Alj^f", 
and F^v = d^A^ - d^A^j, -i[A^, A„]. As usual, the gauge 
invariance in the 5D theory corresponds to the conserva- 
tion of the global symmetry current in the 4D theory. 

At the IR brane we must impose some gauge invariant 
boundary conditions, and we make the simplest choice: 



S = Jd'^x^Tr^^\DX\^ + 3\X\ 



L )zfi 



0. QCD does not a priori fix this 



boundary condition: for example, there may be addi- 
tional terms in the Lagrangian localized at Zm, such as 
and F| + F^. To estimate the sensitivity to such 
terms, we checked that an boundary term with 0(1) 
coefficient (keeping the p mass fixed) yields a 10% cor- 
rection to the p decay constant. We will be using the 
gauge Az ^ and neglecting boundary terms in the La- 
grangian. In this case our boundary conditions are sim- 
ply Neumann. 

The expectation value of the field X is determined by 
the classical solution satisfying the UV boundary condi- 
tion {2/e)X{e) — M for quark mass matrix M: 



Xo{z) 



1 



-Mz + isz^ 



The matrix S is determined by the IR boundary con- 
dition on X. Instead of specifying this condition we 
shall choose E as an input parameter of the model. The 
meaning of E in QCD can be found by calculating the 
variation of the vacuum energy with respect to M |l2l |: 
= {q°'q'^)- We shall assume, as usual, T, — al and 
take M — ruql. 

At this stage the model has four free parameters: m^, 
CT, Zrn and 55. The gauge coupling 55 will be fixed by the 
QCD operator product expansion (OPE) for the product 
of currents, leaving three adjustable parameters. 

We will focus on the Nf — 2 lightest flavors and neglect 
effects of O(m^). Therefore, in Tabled a, /? = 1, 2; a, 5 = 
1,2,3 and t° — cr"/2, where a" are the Pauli matrices. 

Matching the 5D gauge coupling. — We will use the 
holographic duality to relate of the 5D coupling 55 in (0) 
to the number of colors Nc in QCD. The precise sense 
of the holographic correspondence is the equivalence be- 
tween the generating functional of the connected correla- 
tors in the 4D theory Vy4£)[(/)o(x)] and the effective action 
of the 5D theory S5D,cB[(t>ix, e)], with UV boundary val- 
ues of the 5D bulk fields set to the value of the sources 
in 4D theory: 

W4D[0o(a;)] = 5'5D,cff[(/>(a:,e)] at (^(x, e) = (/)o(a;). (4) 

QCD Green's functions can therefore be obtained by dif- 
ferentiating the 5D effective action with respect to the 
sources. In the case that stringy effects can be neglected, 
S5D, cS is simply given by Eq. The action is evalu- 
ated on solutions to the 5D equations of motion subject 
to the condition that the value of each bulk field at the 
boundary 2: = e ^ be given by the source (j) of the 
corresponding 4D operator O (see Tabled). 

We may now fix the 5D gauge coupling by comparing 
the result for the vector current two-point function ob- 
tained from the above prescription with that of QCD. 
Introducing the vector field as V = {Al + A]i)/2, one 
finds, in the Vz{x,z) = gauge, the equation of motion 
for the transverse part of the gauge field: 



dz -dzV^^iq,z) +^V^^iq,z) 



= 0. 



(5) 



Here V^{q,z) is the 4D Fourier transform of V^{x,z). 
The equations of motion are linearized, as is appropriate 
for determination of two-point functions. Evaluating the 
action on the solution leaves only the boundary term 



S 



1 



d^x 



(6) 



(3) 



If V^o^°((?) is the Fourier transform of the source of the 
vector current = q^^t°-q at the boundary then letting 
y(g,2:) = V{q,z)V^{q), we require that F(g, e) = 1. 
Differentiating twice with respect to the source Vb, we 



3 



arrive at the vector current two-point function, 



This allows us to extract the decay constants Fp-. 



nv(-g') 



1 d,V{q,z) 



, (7b) 



where — —q^. For large Euclidean we only need 
to know V{q, z) near the boundary, 



y(g,z) = i + ^in(gV 

which up to contact terms gives 



1 



nv(g') = -2^1nQ 



(8) 



(9) 



On the other hand, we can cornpute Hy from QCD by 
evaluating Feynman diagrams 13<|. The leading-order 
diagram is the quark bubble. 



nv(Q') 



lnQ2 



This leads to the identification 



9l 



(10) 



(11) 



which completes the definition of the action |2Jl . 

Hadrons. — The hadrons of QCD correspond to the nor- 
malizable modes of the 5D fields. These normalizable 
modes satisfy the linearized equation of motion and decay 
sufficiently rapidly near the boundary z — > so as to have 
a finite action. The IR boundary condition gives rise to 
a discrete tower of normalizable modes. The eigenvalue 
of a normalizable mode is the squared mass of the cor- 
responding meson, and the derivative of the mode near 
the UV boundary yields the decay constant. 

To illustrate the above, consider the tower of the p 
mesons. A p wavefunction, ipp{z), is a solution to Eq. 



for an arbitrary component of Vn with q~ 



sub- 



ject to Vp(e) = 0, d;:1pp{Zm) 



and normalized as 



j{dz / z) ij} p{zY — 1. Consider the Green's function cor- 
responding to Eq. (jSJ for an arbitrary component of V^: 



G{q;z,z') 



tppiz)ipp(z') 
p P 



(12) 



(The ie prescription, among other things, guarantees the 
positivity of the spectral function, contrary to the claim 
of Ref. |li|.) One can show that V{q,z') of Eq. (|7b|) is 
given by —{l/z)dzG{q] z, z') at z — e. Now from (|7b|l we 
find: 



nWV) = -;^E 



[V';(e)/e] 



9i V (9^ ~rn^p + ie)ml ' 



(13) 



F'p-iK'ie)/e?^^[^;iO)r, (14) 

95 95 

where Fp is defined by (0| J^lp**) = Fp5°'''ep for a p meson 
with polarization e^. Eqs. © and (|13|) are the holo- 
graphic version of the QCD sum rules. 

In the axial sector (ai and tt mesons), the action to 
quadratic order is 



S ■ 



(Tx 



^9P 



A 



a\2 



(15) 



where we have defined v{z) — niqZ + az'^, A — (Al — 
Ar)/2, and X = exp(i27r"i"). In the Az ^ gauge, 
the resulting equations of motion in 4D momentum space 
are {Ap = Ap± + dpip) 



dz 



-dzK 



fl2 

-^1 



2 2 

95V' 



3 



= 0; 



z 



95V 



+ ^('r'^-^°) = 0; 



-q^dzV" 



9 9 

55 « 



0. 



(16) 
(17) 
(18) 



The oi, being a spin-1 particle, is the solution to Eq. Hlt)|) 
with ■(/'ai(O) — dzipai{zm) = 0. Thc fli decay constant, 
f aj , is given by an expression similar to Eq. H14I) , but 
with p replaced by ai. 

Our theory has all the consequences of chiral symme- 
try built in. Let us derive the Gell-Mann-Oakes-Renner 
(GOR) relation. 



mlf^ = (tou -I- md){qq) = 2mqa. 



(19) 



Since (0|A^|7r) = ifjiq^^ the axial current correlator in 
the = limit has a singularity at = 0: IlA{—q^) 
—fHq^. Using the holographic recipe [cf. Eq. Q], 



2__^ dzA{Q,z) 

Jtt 2 

95 ^ 



(20) 



where ^(0, z) is the solution to Eq. Hlt)|) with q^ = 0, 
satisfying A'{0, z„i) — 0, A{0,e) = 1. The pion is the so- 
lution to Eqs. ((T7|) and (fTH|l . subject to (p'{zm) = fi^) = 
7r(e) — 0. We may construct such a solution perturba- 
tively in ttItt by letting ip{z) — A(0, z) — 1. Then, from 
Eq. H18|l . to leading order in m^. 



7r(z) = ml 



du 



v{uY glu 



-duA{Q,u) 



(21) 



The function /v{u)'^ has a significant support only for 
M ~ Zc = \/'(nqla. The function 9„A(0, u)/{g1u) for such 
small values of u can be replaced by its value at u = e, 
which is related to /tt via H20|l . Performing the integral 



4 



one finds that tt = — mJ/^/(2TOq(T) for z ^ z^- Equations 
(|16|l and H17|l are solved by 1^9 = A(0, z) — \ and tt = const 
for z ^ Zc only if vr = —1, hence mj/^ = 2mgCr + 0(777,^). 

Meson interactions and 5p,rir • — The meson interactions 
can be read from the nonbilinear terms in the 5D action. 
For example, we find that the tt-/? coupling is given by 

9p7rn = 95 J dzipp{z) y—gT^ ^ j ■ (22) 

The normalization of tt is fixed by the pion kinetic term: 
integrating the function in parentheses in Eq. H22|l gives 1. 
One must be aware that this 3-meson amplitude could be 
sensitive to the terms not yet included in our model. 

Predictions. — From Eq. (O and the Dirichlet boundary 
conditions, the p wavefunctions are Bessel functions with 
masses determined by zeroes of Jo[qZm)- Hence, rup = 
2.405/z™ = 776 MeV fixes z„ = 1/(323 MeV). and 
a can then be fit to the experimental values of ttj^ and 
yielding = 2.29 MeV and a = (327 MeYf. These 
parameters correspond to Model A in Table ITU 

1 /2 

The rms error, e^i^s — (X]o('^*^/^)^/") (where 
50 /O is the fractional error of an observable O and n = 4 
equals the number of observables minus the number of 
parameters) for Model A is 15%. 

A global fit to all seven observables (Model B) yields 
the parameters, Zm — 1/(346 MeV), ruq — 2.3 MeV and 
cr = (308 MeV)^ The last column of Table El lists the 
calculated observables in this model. The rms error of 
Model B is a remarkably small 9%. 

Discussion and outlook. — The holographic model of 
QCD studied here is quite crude and depends on only 
three free parameters, but it agrees surprisingly well with 
the seven experimentally measured observables which we 
have studied. There are several ways in which we may 
attempt to extend and improve the model, (i) The glue- 
ball spectrum can be calculated from the gravitational 
and dilaton modes in the theory, which were not in- 
cluded in this study, (ii) It is straightforward to de- 
scribe power corrections in the current correlators |l7l |. 
Here we matched the gauge coupling 35 in our model to 
the leading term — the unit operator — in the OPE of the 
product of currents. Higher dimension operators also ap- 
pear in the OPE, suppressed by powers of the Euclidean 
momentum Q. These corrections can be calculated in 
QCD ||13|. In the holographic model, these corrections 
arise from trilinear and higher terms in the 5D action, 
such as / d^x^X'^F'^. Matching the QCD OPE coef- 
ficients to the coefficients of the 5D action provides a 
method of building and constraining the effective 5D ac- 
tion, (iii) Including the strange quark into the model 
with an approximate SU(3)xSU(3) chiral symmetry is a 
natural extension of the model, (iv) The chiral anomaly 
can be incorporated via a 5D Chern-Simons term, (v) We 
can include corrections to the dimensions of the chiral or- 
der parameters by varying the mass of the corresponding 



TABLE II: Results of the model for QCD observables. Model 
A is a fit of the three model parameters to rriTr, /tt and rUp 
(see asterisks). Model B is a fit to all seven observables. 



Observable 


Measured 
(MeV) 


Model A 
(MeV) 


Model B 
(MeV) 


771 TT 


1 on f?_i_n AAA/i rol 
io9.D±U.UUU4 


1 on t^* 
139.0 


1/11 

141 




775.8±0.5 [8j 


775.8* 


832 




1230±40 [8j 


1363 


1220 




92.4±0.35 [g] 


92.4* 


84.0 


^ P 


345±8 fl5] 


329 


353 


-f ai 


433±13 [6, 16] 


486 


440 




6.03±0.07 [8] 


4.48 


5.29 



fields X in the 5D theory, and we can include running 
of the gauge coupling via logarithmic corrections to the 
AdS geometry. It is interesting to note in this context, 
that those results which follow from partial conservation 
of the axial current, e.g., the COR relation, continue to 
hold as we vary the 5D mass of X in the model 0| . 

This work is supported in part by DOE Grants DE- 
FG0201-ER41195andDE-AC02-76SF00515. J.E. thanks 
the College of William and Mary for support. D.T.S. and 
M.A.S. thank A. P. Sloan Foundation for support. 



[1] J. M. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998). 

[2] A. Karch and E. Katz, JHEP 06 (2002) 043. 

[3] M. Kruczenski et al, JHEP 07 (2003) 049. 

[4] S. Hong, S. Yoon, and M. J. Strassler, JHEP 04 (2004) 

046; hep-th/0409118 hep-th/0410080 
[5] C. Csaki et al, JHEP 01 (1999) 017; L. Girardello et al, 

Nucl. Phys. B 569, 451 (2000); J. Babington et al, Phys. 

Rev. D 69, 066007 (2004); M. Kruczenski et al, JHEP 05 

(2004) 041; H. Boschi-Filho and N. R. F. Braga, JHEP 

05 (2003) 009; G. F. de Teramond and S. J. Brodsky, 

Phys. Rev. Lett. 94, 201601 (2005). 
[6] D. T. Son and M. A. Stephanov, Phys. Rev. D 69, 065020 

(2004). 

[7] M. Piai, A. Pierce, and J. Wacker, |he'p-ph/0405242| 
[8] S. Eidelman et al (PDG), Phys. Lett. B 592, 1 (2004). 
[9] S. S. Gubser, I. R. Klebanov, and A. M. Polyakov, Phys. 
Lett. B 428, 105 (1998). 
[10] E. Witten, Adv. Theor. Math. Phys. 2, 253 (1998). 
[11] I. R. Klebanov and M. J. Strassler, JHEP 08 (2000) 052. 
[12] I. R. Klebanov and E. Witten, Nucl. Phys. B 556, 89 
(1999). 

[13] M. A. Shifman, A. I. Vainshtein, and V. I. Zakharov, 
Nucl. Phys. B 147, 385 (1979). 

[14] N. Mahajan, Phys. Lett. B 623, 119 (2005). 

[15] J. F. Donoghue, E. Golowich, and B. R. Holstein, Dy- 
namics of the Standard Model (Cambridge University 
Press, Cambridge 1992) 

[16] N. Isgur, C. Morningstar, and C. Reader, Phys. Rev. D 
39, 1357 (1989). 

[17] J. Erlich, E. Katz, D. T. Son, and M. A. Stephanov (to 
be published). 



